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Abstract. In this paper I present a new two-dimensional decomposition technique, which models the surface pho- 
tometry of a galaxy with an exponential light profile for both bulge and disk and, when necessary, with a Freeman 
bar. The new technique was tested for systematic errors on both artificial and real data and compared with widely 
used one-dimensional decomposition techniques, where the luminosity profile of the galaxy is used. The comparisons 
indicate that a decomposition of the two-dimensional image of the galaxy with an exponential light profile for both 
bulge and disk yields the most reproducible and representative bulge and disk parameters. 

An extensive error analysis was made to determine the reliability of the model parameters. If the model with an 
exponential bulge profile is a reasonable description of a galaxy, the maximum errors in the derived model parameters 
are of order 20%. The uncertainties in the model parameters will increase, if the exponential bulge function is replaced 
by other often used bulge functions as the de Vaucouleurs law. 

All de compo sition methods were applied to the optical and near-infrared data set presented by de Jong & van der 
Kruit (1994), which comprises 86 galaxies in six passbands. 
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1. Introduction 

The light distribution of disk dominated galaxies is often 
decomposed into a bulge and a disk component which are 
assumed to be physically and dynamically distinct. The 
disk component is flat and governed by rotational dynam- 
ics. The spherical bulge component, though mainly rota- 
tionally supported against gravity, is dynamically a much 
hotter system than the disk. Whether this separation is 
real is hard to say; it is likely that a dynamical interplay 
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exists between the different components in the dense in- 
ner regions of galaxies. Separating both components using 
only the surface photometry of a galaxy has been a long- 
standing problem. 

When we assume that there are distinct components, 
the parameters describing the light distribution of these 
components are of fundamental importance. They reveal 
the common properties among galaxies and especially in 
combination with dynamics and (chemical) content they 
are tracers of galaxy formation and evolution. Two well- 
known relations with fundamental component parameters 
are the constancy of central sur face brightness among 
disks of spiral galaxies (Freeman 197C) and the link be- 
tween bulge-to-disk (B/D) ratio and the Hubble classifi- 
cation sequence. 

Many different decomposition techniques can be found 
in the literat ure (fo r reviews see e.g. Simien 1989, Capac- 
cioli & Caon 1992 ). For most of the methods one postu- 
lates some mathematical functions describing the shape 
of the different components, after which the sum of the 
components functions is fitted to the observed light dis- 
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tribution. Decomposition techniques differ in both the as- 
sumed mathematical functions as well as in the applied 
fitting algorithms. 

The most frequently used function describing the ra- 
dial surface brightness profile of the disks of spiral galaxies 
is the exponential function 



I]disk(r) = Soe-'"/" 
or in magnitudes 
Mdisk(r) = /io + l.086r/h 



(1) 



(2) 



with /io (So) being the central surface brightness (luminos- 
ity) and h the scalelength of the exponential disk. The ex- 
ponential function is not always a good description of the 
disk light profile. After the bulge light is subtracted from 
the luminosity profile, the disk profiles of spiral galax- 
ies sometimes show a deficit of light in their inner regions 
(Type II profiles, Freeman 1970). An improved disk model 
in the form of a modified exponential profile (Kormendy 



1977) is therefore sometimes fitted to these Type II pro- 



files. Furthermore, it should be noted that the exponen- 
tial model profile might need a truncation. Especially in 
edge-on galaxies a sudden decrease in light has often been 
detected at large radii (van der Kruit 1979 ). 

Compared to the disk profile there is less consensus 
on the mathematical function to be used for the bulge 
light profile, because only the central bright part of the 
bulge can be seen in most spiral galaxies. Away from the 
center the bulge light is hidden underneath the disk light 
for a face-on galaxy and therefore only the central region 
can be used to determine the shape of the bulge lumi- 
nosity profile. The most widely used bulge function is the 
r^/'* law (de Vaucouleurs 1948). This function generally 
gives a good description of the light distribution of el- 
liptical galaxies and was first used for spiral bulges, in 
combination with an exponential disk, by de Vaucouleurs 
( |1959| ). The motivation to use the same fitting function 
for elliptical galaxies and bulges is an assumed evolution- 
ary or at least structural sequence from ellipticals to spi- 
rals (or vice versa) . The shape of elliptical luminosity pro- 
files are not completely undisputed and therefore all light 
profil e func tions proposed for ell ipticals (e.g. Hu bble p ro- 
files ( |1930| ), King profiles (|1966D , Jaff e pro files ( |l983|) or 
ge neraliz ed exponentials (Caon et al. 1993| ; D'Onofrio et 
al. 1994)) cou ld in principle also be used for bulges. Kor- 
mendy (1977) showed that the parameters of the Hubble, 
King and de Vaucouleurs profiles describe the same phys- 
ical quantities for ellipticals. 

One should be careful in making the link between ellip- 
tical galaxies and bulges of spiral galaxies. The presence of 
the disk will influence the dynamics of the bulge, especially 
for late- type systems with a low B/D ratio. Bulges are 
for a considerable fraction rotationally supported against 
gravity and not largely pressure supported as are elliptical 
galaxies (Kormendy & lUingworth 1982; Kormendy 1993 



and references therein). Also deviations from the r^/^ law 
have bee n obse rved in bulges of edge- on spi rals (Fran kston 
& Schild |l976i K ormen dy & Bruzual |l978|, Bu rstein |l979 , 
Jense n fc T huan |l979i Shaw & Gilmore |l989| , Wainscoat 



et al. 1989 Kent et al. 1991). These deviations have mo- 



tivated Frankston & Schild (^976|) , Kent et al. ( |l99l| ) and 



Andredakis & Sanders (1994) to propose exponential func 
tions for bulge profiles. 

Decomposition techniques using the change in ellip- 
ticities and position angles of the isophotes (Kent p_986 ) 
have the advantage that they do not have to assume fitting 
functions. They can work perfectly well assuming that spi- 
ral galaxies are only made of ellipsoids/tori with changing 
inclinations. But as real galaxies contain bars, spiral arms 
and dust lanes these methods will have systematic errors. 
They only work reasonably well on systems with a high 
inclination, where the difference in the flattening of the 
disk and the bulge component is easily measurable. 

Decompositions using all the pixels in the full image 
instead of the one-dimensional (ID) profile of galaxies take 
little precedence in the literature. The main reason for 
this being the lack of computer power and an additional 
reason being the difficulties due to the presence of bars 
and spiral arms. These are conveniently averaged out in 
ID profiles. Two-dimensional (2D) fitting techniques have 
been applied before on small samples of elliptical galaxies 



and SO's (Capacc ioli et al. |1987t Simien & Michard |1990 ; 
Scorza & Bender 1990| ), b ut on ly once before on a large 
set of spiral galaxies (Byun |1992 ). Two-dimension al fitt ing 
has the same advantage as Kent's method (Kent |1986 ) in 
that one uses the difference in projected ellipticities of 
disk and bulge. The sample used here was selected to be 
face-on and the differences in ellipticities are expected to 
be small. Still the 2D fitting technique is applied, because 
it has the advantage to ID fitting that non-axisymmetric 
components can be fitted as well. The method applied 
here has a non-axisymmetric component in the form of a 
bar, which will improve the fitting results of the disk and 
especially the bulge component. 

The bulge and disk parameters determined by the 2D 
method will be used in subsequent papers to determine the 
relationships among the structural parameters of galaxies. 
To assess the reliability of these relationships a thorough 
error analysis is needed. In this paper I will discuss several 
sources of error, most notably the effects of ID versus 
2D fitting, the uncertainty in the shape of disk and bulge 
profiles, influences of measurement errors and the effect of 
different radial weighting functions. 

The structure of this paper is as follows. The observa- 
tional data are briefly described in Sect. |. The 2D fitting 
technique is explained and tested in Sect. |. In Sect. ||the 
2D fitting technique is compared with several ID fitting 
techniques. The results of all tests and comparisons are 
used in the error discussion in Sect. |[ The main conclu- 
sions are summarized in Sect. O- 
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2. The data 

In order to examine the parameters describing the global 
structure of spiral galaxies, 86 systems were observed in 
the B, V, R, I, H and K passbands. A full description of 
the observations and data reduction can be found in de 
Jong & van der Kruit ( 1994 , hereafter Paper I) , which 



will be repeated only very briefly here. The galaxies in 
this statistically complete sample of undisturbed spirals 



were selected from the UGC (Nilson 1973) to have red 
diameters of at least two minutes of arc and minor over 
major axis ratios larger than 0.625. The galaxies were im- 
aged along the major axis with the Im Jacobus Kapteyn 
Telescope on La Palma in the B,V,R and / passbands and 
with the United Kingdom Infra-Red Telescope on Hawaii 
in the H and K passbands. Standard reduction techniques 
(bias subtraction, flatfielding by twilight flatfields, calibra- 
tion by standard stars) were used to produce calibrated 
images. The sky brightness was determined outside the 
galaxy in areas free of stars and its uncertainty (mainly 
due to flatfield limitations) constitutes one of the main 
sources of error in the derived parameters. 

The ellipticity and position angle (PA) of each galaxy 
were determined at an outer isophote. The radial surface 
brightness profiles were determined by calculating the av- 
erage surface brightness on elliptical rings of increasing 
radii using the determined ellipticity and PA. Internal and 
external comparisons showed that the derived parameters 
were well within the estimated errors. These estimated 
errors are included in the analysis discussed here. 

3. Two-dimensional decomposition 

In this section the 2D fitting technique is described. The 
motivation for using the 2D method was the large number 
of galaxies in the sample with a pronounced bar, which 
can not be fitted in ID models. The different model com- 
ponents are described and the fitting procedure followed is 
explained in some detail. The fitting technique was tested 
on both artificial and real data and is shown to be very 
accurate in most realistic cases. Finally, the results for the 
data set are presented for the B and the K passband. 

3.1. Advantages of two-dimensional fitting 

In the literature one encounters mainly the use of ID 
profiles to perform bulge/disk decompositions. The ex- 
traction of the profiles improves the signal-to-noisc, but 
the non-axisymmetric information present in the image is 
lost. Bars can have considerable influence on proflles (see 
for example in Paper I UGC 89, UGC 6536, UGC 7523, 
UGC 7594, UGC 8865 and UGC12776). These features 
will make ID bulge/disk decompositions incorrect, even 
if the flts seem correct and the values are low. One- 
dimensional models fitted to azimuthally averaged profiles 
can never include a bar. 



A considerable fraction of spiral galaxies are barred. 
Of the 86 galaxies of our sample only 13 were classi- 
fic d as n on-barred according the RC3 (de Vaucoulcurs et 
al. 1991 ) and 12 galaxies had no bar classification. There- 
fore fitting bars is desirable, especially if one considers 
that bars are more pronounced in the near-IR (Block & 
Wainscoat 1991 ). This can only be done by fitting the full 
2D images. An additional advantage of the 2D technique 
is that the difference in the flattening of the bulge and the 
disk is also used for the few systems in the sample with 
higher inclination. 

The technique of fitting models to the full images of 
spiral galaxies has little precedence in the literature. Shaw 
& Gilmore ( |1989| ) fitted 2D models to two edge-on galax- 
ies, the configuration where disk and bulge are the most 
distinct. Byun ( |1992| ) used a data set of 1355 / passband 
images of galaxies with inclinations larger than i > 40°. 
From tests on artificial data, Byun showed that the 2D 
method was better in reproducing the model input pa- 
rameters than the ID method. Neither Shaw & Gilmore 
nor Byun included a bar in their fits. 

3.2. The model components 

The 2D model consists of two or three components. First 
of all, a spherically symmetrical bulge with an exponen- 
tial radial light distribution. The use of an exponential 
bulge profile instead of the more widely used r^/* law pro- 
file (de Vaucoulcurs 1948) wa s motivated by the work of 
Andredakis & Sanders (1994) and the work presented in 
Sect. ^. Bulge parameters are normally expressed in effec- 
tive parameters which translates the exponential law into 



I]buigc(r)=See-i-6^9('-/--i), 



(3) 



where the effective radius (ro) encloses half the total lu- 
minosity and Se is the surface brightness (/ie in mag) at 
this radius. In the cartesian coordinate system of the CCD 
image r should be read as \/ x'^ +2/^, with the center of 
the coordinate system at the galaxy center. 

The second component, the disk, is described 
by Eq. (|^) and has the usual two free parame- 
ters of an exponential light distribution (Eq and h). 
Due to inclination though, the disk has two addi- 
tional parameters: minor over major axis ratio (b/a) 
and position angle (PA) and r should be read as 



([a:cos(PA) 2/sin(PA)]f )2 -|- (-xsin(PA) 2/cos(PA))2 
in the carthesian coordinate system. Fitting with b/a and 
PA as free parameters was tried, but it turned out that 
the fitting routine often adjusted the b/a and PA in such 
a way that spiral arms and bars were modeled, instead 
of the global disk properties. Therefore I decided to keep 
b/a and PA fixed to the values determined at the outer 
isophotes (the same values that were used for extracting 
the radial profiles). 
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When the galaxy image showed a clear bar and when 
this bar resulted in an identifiable feature in the luminos- 
ity profile, a bar was added as a third component to the 
model. A Freeman bar ( 1966 ) was used, which is one of the 
few available analytic 2D descriptions for the luminosity 
of a bar 



(4) 



where the free parameters are So, bar, the bar central sur- 
face brightness, and Obar and &bar, the semi major and mi- 
nor axis of the bar respectively. Such a bar has elliptical 
isophotes and has its position angle (PAbar) as additional 
free parameter. No inclination dependent corrections were 
applied to the bar profiles, since the studied galaxies are 
not inclined very much. 

Observations of galaxies are distorted by seeing. The 
model light distributions have to be corrected for this ef- 
fect. Seeing is only important at the center of the galaxy 
where the light distribution strongly peaks. The disk and 
bar light distributions are far less peaked and generally do 
not dominate the center compared to the bulge, therefore 
only the bulge model profile was corrected for seeing. 

To account for the seeing effects, the model profiles of 
the bulge were convolved with a Gaussian Point Spread 
Function (PSF) with the same dispersion (cr) as the Gaus- 
sians that were fitted to some field stars in the frame. For 
a radially symmetric light distribution around the center 
the seeing convolved profiles are described by 



T.{x)lQ{xr/(T'^)e''^/'^''\dx, (5) 



where is the intrinsic surface brightness profile, cr 

the dispersion of the Gaussian PSF and /o the zero-order 
modi fied B essel function of the first kind (Pritchet & 
Kline [19811 ). 

The total model is simply the sum of the three com- 
ponents 

Stot(a;, y) = Sdisk(a;, y) + Ebuigo(a;, y) + Sbar(a;, y). (6) 
3. 3. Fitting procedure 

To fit models to the data points a non-linear fitting al- 
gorithm capable of accepting different weights for each 
data point was applied. Non-linear fitting algorithms are 
particularly sensitive to the initial values provided, when 
searching for the minimum in the reduced of the fit. If 
the initial values are not "reasonable" , the fitting program 
can end up in a wrong local minimum. The results of the 
ID decompositions described in Sect. ^ were used as ini- 
tial values. The initial values for the bar were estimated by 
eye. The routine generally converged to the same result to 
within the formal errors, independent of the initial values. 
The formal fit errors were usually much smaller than the 
errors due to the uncertainties in the measurements and 
these errors will be discussed in detail in Sect. ^ and ||. 



In decomposing ID profiles it is common practice to fit 
in the logarithmic (magnitude) regime. This means that 
one is effectively trying to minimize the relative errors be- 
tween model and data. Minimization in the logarithmic 
regime is not possible in the 2D case; because of the noise 
some pixels will have negative values (below sky level) in 
the outer parts of the galaxy image. To minimize relative 
errors in the linear regime, a difference between model and 
data in a low surface brightness region has to be given 
much more weight during the fitting than the same differ- 
ence in a high surface brightness area. A weight function 
of the form e'"/'* was used, where h is the initial estimate 
of the disk scalelength and r is the inclination corrected 
distance of the pixel from the center. To reduce comput- 
ing time all pixels outside 2.5 initial disk scalelengths were 
averaged over 5x5 pixels and given proportionally more 
weight. 

The 2D fitting procedure consisted of several steps. 
First, all images of a galaxy obtained in different pass- 
bands were aligned, freed of foreground stars and their 
center was determined from the R passband image (see 
Paper I). This center was fixed while fitting the model 
components to the data. After the fitting routine had con- 
verged, the model light distribution was subtracted from 
the data. The points in the difference image that deviated 
more than 6 sigma (missed cosmic ray events and faint 
stars) were flagged and not used in the next iteration. 
This process was repeated twice, each time taking the ini- 
tial estimates and the scalelength of the weight function 
from the results of the previous step. In general, the rou- 
tine had already converged to a satisfactory result after 
the second step. 

3.4- Tests on artificial data 

The 2D fitting routine was extensively tested on artificial 
images to determine its reliability. Artificial images are 
not really representative of true galaxies, but can give an 
indication of the systematic effects due to measurement 
errors. By varying one by one the observables in the arti- 
ficial images one can investigate which measurement error 
influences a model parameter the most. 

The artiflcial images had the characteristics of a typ- 
ical R passband observation: all artiflcial galaxies had an 
exponential disk with a //q equal to 20 mag arcsec"^, a 
scalelength of 20" and axial ratio of 0.75. Just as with 
the typical observation the pixel size was set at 0.3", the 
seeing at 1.5" FWHM and the sky surface brightness at 
20 mag arcsec"^ in the artificial images. An exponential 
bulge was added to each image, with the bulge param- 
eters chosen from a range in effective surface brightness 
and radius . The images were created with photon and 
read out noise and a few areas were set to undefined val- 
ues to mimic the removal of foreground stars. The initial 
fit estimates for the disk and bulge parameters were set 
at 10-30% off the intrinsic model values to check conver- 



R.S. de Jong: Near-IR and optical observations of 86 spirals. II Determination of bulge and disk parameters 



5 



Table 1. The relative errors between input artificial image parameters and fitted bulge and disk parameters determined using 
the 2D fit method. Tabulated are the results for the standard artificial galaxy as described in the text, as well as the effect of a 

wrong estimate of one of the parameters which were normally kept fixed to the observed values while fitting the model (see notes). 
Relative errors are listed, with A/io = Mo, art — Mo, fit and A/ie = Pe,art — Me,flt in mag arcsec"^. Ah = 2(/iart — /iflt)/(/iart + ftflt) 
and Are = 2(re,art - ?'e,fit)/(^e,art + r-e,fit) dimensionless. 
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Notes: 



(1) The bulge effective surface brightness of the artificial galaxy in mag arcsec"^. 

(2) The bulge effective radius of the artificial galgixy in arcsec. 

(3) Fit results using the standard artificial image (see text) with //o = 20 R-mag arcsec"^ and h = 20". 

(4) As (3), but artificial image and fitting model have seeing of 2.5" instead of 1.5". 

(5) As (3), but fitting model has 1.4" seeing instead of the 1.5" of the artificial image. 

(6) As (3), but fitting model has 1.6" seeing instead of 1.5". 

(7) As (3), but with sky level of the fitting model 1% too low compared to the real value in the artificial image. 

(8) As (3), but with sky level of the fitting model 1% too high compared to the real value in the artificial image. 

(9) As (3), but with b/a = 0.65 in fitting model and 0.75 in the artificial image. 

(10) As (3), but with b/a = 0.85 in fitting model and 0.75 in the artificial image. 

(11) As (3), with error of 20° in PA between image and model. 
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Table 2. The test results of the 2D fit using an R passband image of UGC 438. Before the test a bulge with /^o=18.954 i?-mag 
arcsec"'^ and re=2.5" (the parameters resulting from the initial fit) was subtracted from the original image, leaving a bulgeless 
test image. The removed bulge was for the test replaced by an exponential bulge of the indicated model parameters (model /ic, 
model re) and fitted by the 2D fit routine. The resulting fitted /io and /Xo are in _R-mag arcsec"'^ and the h and rc are in arcsec. 
NC indicates that no fit convergence was reached. 
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gence. The final results showed very little sensitivity to 
the initial estimates. 

The effect of different bulge-to-disk ratios was tested 
using the artificial images with the different bulge parame- 
ter values. The relative differences between artificial image 
input parameters and the model fit output parameters are 
listed in Table 0. The fitting routine reproduced the bulge 
and disk parameters to a very high degree of accuracy. 
The disk parameters were reproduced with less than 1.5% 
error. The bulge parameters were less well reproduced, 
especially for very small effective radii and low effective 
surface brightnesses. These bulges are so small that they 
are unresolved with 1.5" seeing in the small region where 
they are dominating over the disk surface brightness. For 
these extreme cases it is better to convolve not only the 
bulge model with the seeing, but also the disk. Increasing 
the seeing to 2.5" in both the artificial image and the fit- 
ting model had no effect on the determination of the disk 
parameters and slightly decreased the fit performance for 
the bulge parameters, as can be seen in Table |l[ 

Some observables measured from the images (seeing, 
sky background level, b/a and PA) were used as fixed pa- 
rameters in the fit models. To estimate their relative ef- 
fects on the derived parameters, all artificial images were 
fitted again, but the fixed parameters in the fitting model 
were given a wrong value with respect to the true values in 
the artificial image. The typical maximum error was used 
for each of the fixed parameters. The test results can be 
found in Table |l|. From these tests the largest error in /io 
is expected to result from a wrong estimation of the sky 
background. The error will be of order 0.1 mag arcsec"^, 
increasing, of course, for lower surface brightness galaxies. 
The errors in the scalelength are dominated by errors in 
sky brightness and 6/a, and both result in errors less than 
10% in the typical case. The errors in the bulge param- 
eters are dominated by their brightness and scalelength 
relative to the disk parameters. The parameters of the 
brighter bulges are most effected by errors in the seeing. 
The effective radii of the lower surface brightness bulges 
are also influenced by sky brightness errors. Except for the 



extreme low surface brightness bulges and for bulges with 
effective radii smaller than 1" the errors are never larger 
than 20%. 

The fitting routine was also tested on some artificial 
images with a bar, but with eight free parameters only 
a limited parameter space could be explored. In general, 
the routine converged to satisfactory results in just a few 
iterations. 

3.5. Tests on UGC 438 

Clearly, artificial images are not a good representation of 
real galaxies. Therefore, the fitting routine was also tested 
for different bulge-to-disk ratios on an R passband im- 
age of UGC 438. A 2D fitted bulge was subtracted from 
the original image and replaced by artificial exponential 
bulges with a range in effective surface brightness and ra- 
dius. The fit results for the different bulge parameters are 
listed in Tabic |. 

The fitting routine generally reproduces the parame- 
ters quite well. The disk parameters change little for the 
different bulges (/io changes at most by 0.04 mag arcsec"^ 
and h by about 4%), but as with the tests on artificial 
images, the fitting routine breaks down for bulges with a 
small effective radius and a low surface brightness. Next 
to the explanations for this effect mentioned earlier for the 
artificial images, one has another problem when one tests 
on real images. The model bulge that was originally sub- 
tracted to create a bulgeless image might not be perfectly 
correct. One is then left with some residuals which have 
the greatest infiuence when later replaced by a low surface 
brightness bulge. 

3. 6. Resulting parameters 

The results from the 2D fits for all galaxies in our sample 
are presented in Table || for the B and the K passband. 
This table contains observed data values. No corrections 
were applied for inclination and for internal and Galactic 
extinction. In 23 out of 86 cases it was found that the fits 
of at least the bulge improved if a bar was also included. 
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Table 3. The results from the 2D fit model with fxo, He and /ito.bar in mag arcsec ^, h, Ve, Obar and 6bar in arcsec and PAbar in 
degrees. NP means a non-photometric observation, but the scale parameters are still determined. 
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24.0 


20.47 


2.1 


22.15 


13.2 


4.9 


19.6 


17.64 


16.2 


15.67 


2.0 


18.06 


16.5 


6.2 


22.5 


1551 


22.47 


25.8 


24.81 


2.6 










18.98 


24.4 


20.06 


2.9 










1559 


22.48 


20.0 


23.42 


3.8 










20.08 


19.6 


20.21 


5.8 










1577 


22.44 


22.1 


20.89 


2.2 


22.80 


32.9 


10.8 


144.5 


18.26 


17.6 


16.01 


1.8 


18.42 


26.8 


8.9 


145.1 


1719 


22.45 


21.7 


20.99 


1.6 


22.88 


13.3 


5.2 


127.6 


17.73 


14.3 


16.68 


1.9 


18.45 


15.9 


5.2 


129.3 


1792 


21.65 


17.1 


21.31 


1.3 










17.48 


14.4 


17.21 


1.8 










2064 


22.28 


20.1 


20.78 


0.8 










18.01 


17.6 


16.25 


1.1 










2081 


22.31 


19.1 


23.46 


1.4 










19.44 


19.4 


19.83 


2.5 










2124 


22.34 


23.0 


20.07 


2.7 


21.96 


21.4 


7.0 


100.1 


18.11 


21.3 


15.89 


2.7 


17.62 


21.3 


6.1 


100.2 


2125 


23.20 


26.8 


21.01 


1.7 


23.47 


34.6 


6.9 


79.5 


















2197 


22.57 


18.2 


23.43 


2.1 










17.97 


12.3 


18.60 


2.1 










2368 


23.67 


28.3 


19.82 


1.3 


22.50 


52.9 


14.0 


156.9 


18.88 


18.1 


15.19 


1.7 


18.41 


48.0 


11.9 


156.9 


2595 


































3066 


22.03 


14.6 


21.40 


0.4 










17.09 


11.2 


17.25 


1.6 










3080 


21.99 


17.2 


19.88 


0.2 










18.21 


15.1 


18.79 


1.7 










3140 


20.90 


13.1 


20.46 


2.0 










16.96 


11.3 


15.99 


2.1 










4126 


21.83 


22.6 


20.01 


1.4 


23.34 


45.7 


9.6 


38.1 


18.02 


18.7 


16.09 


1.8 


19.01 


35.2 


10.8 


48.5 


4256 


21.18 


17.6 


19.22 


0.7 










17.27 


15.8 


16.15 


2.1 










4308 


21.34 


20.0 


19.78 


0.8 










17.66 


16.4 


17.06 


1.9 










4368 


21.52 


17.6 


22.52 


2.0 










17.81 


14.7 


18.89 


2.2 










4375 


21.31 


20.2 


22.18 


2.1 










17.21 


16.9 


18.46 


2.5 










4422 


22.04 


31.0 


20.92 


3.5 










18.38 


27.0 


16.55 


3.4 










4458 


21.72 


19.8 


19.50 


3.0 










17.75 


17.9 


15.05 


2.5 










5103 


20.50 


17.8 


19.84 


1.3 










16.56 


14.5 


16.07 


1.9 










5303 


21.32 


36.8 


21.41 


3.0 










17.89 


37.7 


17.72 


4.3 










5510 


20.66 


20.4 


19.22 


1.3 










17.37 


18.4 


16.53 


2.2 










5554 


20.98 


19.1 


19.31 


1.5 










16.94 


16.7 


15.66 


2.1 










5633 


23.13 


27.8 


23.67 


0.9 


22.99 


23.5 


3.0 


161.5 


20.01 


29.4 


21.03 


1.5 


19.88 


23.2 


3.7 


162.5 


5842 


21.44 


30.3 


21.01 


0.9 










18.01 


27.4 


18.52 


2.2 










6028 


20.49 


15.3 


21.57 


0.5 










17.27 


13.0 


18.66 


1.6 










6077 


20.88 


18.4 


20.17 


1.1 










17.49 


18.4 


17.22 


2.1 










6123 


NP 


27.1 


NP 


2.1 










17.29 


25.0 


15.46 


2.0 










6277 


20.88 


27.9 


18.98 


1.3 


— 


— 


— 


— 


17.25 


21.5 


16.54 


2.9 


— 


— 


— 


— 


6445 


20.55 


16.1 


20.19 


2.6 


20.94 


11.6 


4.1 


157.6 


16.83 


14.0 


16.17 


2.6 


17.11 


11.4 


4.5 


154.4 


6453 


20.91 


21.3 


22.04 


7.9 










17.16 


16.3 


17.87 


5.6 










6460 


20.66 


27.6 


19.63 


1.8 










17.15 


24.1 


16.53 


2.1 










6536 


NP 


21.7 


NP 


1.9 


NP 


13.5 


9.3 


9.0 


18.57 


20.3 


15.47 


1.9 


17.80 


12.7 


8.1 


16.5 


6693 


21.60 


18.9 


21.72 


0.7 










17.99 


16.7 


17.13 


1.3 










6746 


21.38 


17.7 


20.18 


2.2 










17.30 


16.0 


15.64 


2.1 










6754 


22.71 


32.0 


21.58 


4.2 










18.32 


21.7 


15.77 


1.8 










7169 


20.11 


13.2 


20.20 


1.7 










16.55 


10.0 


16.38 


2.3 










7315 


19.99 


14.8 


20.75 


1.5 










16.03 


14.6 


16.96 


2.1 










7450 


21.17 


64.9 


20.14 


7.8 










17.30 


48.5 


15.72 


6.3 










7523 


21.49 


36.0 


19.61 


4.0 


22.10 


47.2 


11.8 


144.5 


17.70 


32.5 


15.11 


3.1 


17.87 


44.1 


12.9 


143.6 


7594 


21.03 


52.4 


19.16 


3.6 


21.72 


50.2 


21.5 


11.8 


16.94 


44.8 


15.21 


4.1 


17.73 


51.2 


19.6 


8.7 
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Table N. -continued. 



UGC 



B 



K 





MO 


h 






/iO,bar 


f^bar 


6bar 


PAbar 


Mo 


h 




r-e 


Mo, bar 




6bar 


PAbar 


7876 


NP 


21.2 


NP 


1.7 














18.55 


22.1 


20.06 


4.4 














7901 


20.07 


25.1 


20.06 


3.0 














16.06 


19.7 


16.22 


3.8 














8279 


20.52 


13.6 


22.87 


7.0 














16.90 


11.5 


18.59 


5.9 














8289 


21.80 


33.2 


19.59 


3.2 














17.90 


21.1 


15.84 


3.2 














8865 


21.89 


32.0 


20.10 


2.7 


21.83 


19.8 


9.4 


169.2 


18.27 


29.3 


15.65 


2.3 


17.84 


27.2 


9.5 


168.4 


9024 


24.08 


28.2 


21.64 


2.1 


— 


— 


— 


— 


22.07 


45.4 


19.15 


3.6 


— 


— 


— 


— 


9061 


22.63 


62.0 


20.08 


2.1 


21.97 


15.0 


6.5 


111.3 


18.75 


39.3 


16.22 


2.8 


17.73 


16.8 


4.0 


110.2 


9481 


21.22 


17.6 


21.36 


1.1 


22.87 


27.0 


4.3 


94.2 


18.03 


17.5 


17.26 


1.6 


18.46 


21.5 


6.2 


98.7 


9915 


NP 


17.7 


NP 


2.6 


— 


— 


— 


— 


17.21 


15.7 


16.58 


2.3 


— 


— 


— 


— 


9926 


20.13 


17.9 


20.16 


3.1 


— 


— 


— 


— 


16.45 


16.6 


15.96 


4.3 


— 


— 


— 


— 


9943 


20.40 


20.0 


20.92 


2.7 


21.37 


19.8 


6.2 


73.3 


16.60 


17.6 


16.55 


2.3 


17.16 


19.7 


5.9 


70.3 


10083 


21.51 


23.3 


22.65 


4.3 


22.08 


44.8 


8.8 


149.2 


17.60 


18.9 


18.64 


4.2 


18.50 


44.2 


10.1 


149.3 


10437 


24.03 


33.6 


23.28 


8.9 


— 


— 


— 


— 


NP 


17.0 


NP 


7.1 


— 


— 


— 


— 


10445 


21.76 


20.3 


22.99 


2.3 


— 


— 


— 


— 


19.00 


19.6 


20.08 


4.5 


— 


— 


— 


— 


10584 


21.77 


23.8 


21.22 


1.7 


— 


— 


— 


— 


18.18 


19.4 


17.47 


2.2 


— 


— 


— 


— 


11628 


22.27 


38.8 


20.10 


3.1 


— 


— 


— 


— 


17.19 


21.9 


15.09 


2.6 


— 


— 


— 


— 


11708 


21.51 


17.9 


21.77 


1.9 


— 


— 


— 


— 


NP 


14.4 


NP 


2.0 


— 


— 


— 


— 


11872 


20.46 


18.3 


20.68 


7.0 


20.61 


13.2 


3.7 


17.6 


15.71 


12.6 


14.94 


2.5 


16.01 


12.8 


6.0 


13.3 


12151 


23.27 


26.1 


24.78 


6.4 


— 


— 


— 


— 


20.14 


21.7 


21.19 


4.8 


— 


— 


— 


— 


12343 


21.95 


50.7 


21.32 


4.3 


22.01 


66.9 


11.0 


7.0 


17.67 


39.7 


15.94 


3.4 


17.67 


68.2 


12.6 


10.7 


12379 


21.99 


19.3 


20.40 


2.1 










17.44 


15.7 


15.33 


1.9 










12391 


21.50 


15.7 


18.03 


0.2 










17.79 


13.8 


17.36 


1.2 










12511 


22.46 


24.4 


22.38 


1.6 










18.35 


13.2 


18.09 


2.1 










12614 


20.94 


21.6 


18.76 


0.7 










17.24 


20.2 


15.00 


1.2 










12638 


22.17 


21.4 


21.74 


1.3 










18.27 


19.8 


18.24 


2.6 










12654 


21.76 


19.6 


22.73 


2.0 










17.98 


15.2 


18.82 


3.0 










12732 


23.66 


29.8 


25.30 


11.0 










NP 


37.5 


NP 


10.5 










12754 


21.81 


52.7 


22.21 


5.2 


21.83 


48.1 


6.7 


98.0 


NP 


49.2 


NP 


7.6 


NP 


44.6 


8.3 


100.6 


12776 


23.39 


61.3 


20.01 


2.0 


22.02 


22.6 


8.9 


171.5 


19.30 


30.2 


15.52 


2.0 


17.90 


21.0 


7.9 


171.6 


12808 


20.31 


13.8 


17.50 


0.7 


























12845 


22.69 


24.6 


23.32 


2.4 










NP 


18.1 


NP 


2.5 











Figure |i| shows an example of an original, a model and 
a model subtracted image of UGC 89. The relative dif- 
ferences between model and data are quite small and the 
resulting difference image can be used to study small scale 
structures in the bulge and bar region and to study spiral 
arms. UGC 89 is a good example of a galaxy that can only 
be fitted correctly by a 2D method with a bar, as can be 
seen more clearly in the resulting ID profile of Fig. ^. 

In testing both artificial and real images it appears 
that the 2D fitting routine has difficulties fitting bulges 
with a low surface brightness and a small effective radius 
compared to the disk parameters. There are few galaxies 
in the sample with a fitted < 1. This is only found 
in the B and V passbands when galaxies have Type II 
profiles. The routine fits a small bulge to these galaxy 
images to avoid filling up the "central hole" in the disk. 
The Type II behavior is always less pronounced at longer 
wavelengths and indeed in the K passband there are no 
galaxies with ro < 1. Fortunately, there are no galaxies in 
the sample with a small effective radius and a low effective 
surface brightness of the bulge relative to the disk central 



surface brightness. Using these limitations in parameter 
space, the standard column in Table |l| indicates that the 
errors intrinsic to the fitting routine are expected to be 
less than 10%. The errors caused by other sources will be 
discussed in the next sections. 

4. Comparison of different decomposition methods 

The structural parameters derived in the previous section 
will be used in subsequent papers of this series. In order 
to assess the reliability of the derived relations in these 
papers, a good estimate of systematic and random errors 
in the fitted parameters has to be determined. The 2D fit- 
ting results will be compared to various conventional ID 
methods to demonstrate the increased accuracy. The ID 
methods will also be used to investigate the two most im- 
portant sources of error in the bulge and disk parameters, 
namely the sky background error and the uncertainty in 
the shape of the bulge profile. The 2D method was too 
time consuming to be used for these tests. 
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Fig. 1. Left: UGC 89 R passband image, with isophotes overlaid from 18 to 24 _R-mag arcsec"^ in steps of 0.5 mag arcsec"^. 
Center: Fitted model on the same grayscale and contour level as the left image. Right: Residual image using the 2D fit model 
of the center image. The ellipse indicates the area used for the fit and the b/a and PA of the fitted exponential disk. The fitted 
bar component shows up as the inner ellipse in the residuals. Structure in the bar region and the two arms coming off the ends 
of the bar are clearly visible and can not be fitted with this simple model. 



4-.1. One- dimensional decompositions 

One-dimensional decomposition methods are well 



known in the literat ure ( Kormendy 1977; Schombert & 



Bothun |1987| ; Simie n |l989| ; Capaccioli & Caon [l992| ; An 



dredakis & Sanders 1994) and will be described briefly. 
The resulting parameters are only available in electronic 
form. 

The most elementary way to obtain the disk param- 
eters, the "mar king the disk method" , was also used by 
Freeman ( 197C ), when he found the disk central surface 
brightness to be constant among galaxies. The linear part 
of the luminosity profile, plotted on a magnitude scale, 
was marked and a linear least squares fit was made to the 
data points in the indicated range. To be able to compare 
the disk parameters of a galaxy in different passbands, 
I used the same range in radii for all passbands. The 
resulting parameters can be quite sensitive to the mini- 
mum and maximum radius chosen to fit. The difficulties 
of the "marking the disk" method have been discussed 
by Giovanelli et al. (1994) and results obtained by differ- 



ent authors were compared by Knapen & van der Kruit 
(1991), showing remarkable differences. These differences 
were mainly caused by a change in scalelength at a lower 
surface brightness. 

The "marking the disk" method yields disk parame- 
ters that are intuitively correct for the human eye. How- 
ever, the luminosity profile is a combination of bulge and 
disk light and to get correct results both should be fitted 



simultaneously (Kormendy 1977). The numerical method 
to decompose the luminosity profiles used here is essen- 
tially identical to the method described by Andredakis & 
Sanders (1994). A non- linear minimalization routine 



was used to fit the model profiles to the data points in the 
logarithmic regime. Both bulge and disk model profiles 
were convolved with the seeing PSF using Eq. (||) and the 
fits were limited by the same maximum radii used for the 
"marking the disk" fit. 

As already indicated in Sect. |^, several different ra- 
dial luminosity laws have been proposed for the light of 
bulges. One of the more general forms for the light profile 
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of the bulge is the genera lized exponential law, originally 
proposed by Sersic (1968): 



E(r) = Soe 



-{r/hy 



(7) 



This generalized exponential law has been applied to fit 
elliptical and SO galaxies by Caon et al. ( 1993D and to 
dwarf ellipticals by Young & Currie (1994). 

As a first step, the fitting of an exponential disk and 
a generalized exponential bulge profile was tried. A wide 
variety of initial values for the five free parameters were 
tried, but for most of the galaxies the fits did not con- 
verge to physically acceptable values (negative Sq, Sc, h 
or re). The value of n in Eq. (Q) was of order 0.5-5 in the 
cases where the fit did converge, that is for the galaxies 
with pronounced bulges. For most galaxies the bulge light 
dominates over the disk light at only a few data points 
and these few points do not carry enough information to 
limit the shape parameter n. 

To reach more stable results the same fits were made 
again with n fixed to values 1, 2 and 4. With n = 4 and af- 
ter redefining Eq and h into effective parameters, Eq. (^ 
translates in the most commonly used bulge fitting func- 
tion, the de Vaucouleurs (1948) or r^^'^ law: 



E(r) = See- 



-7.67(r/r„i''*-l) 



(8) 



Setting n = 2 in Eq. (|^) gives an "r^/^ law" profile: 

S(r) - Eee-^-^''(''/''='''-i^ (9) 

In the case of rt = 1 in Eq. one has the exponential 
law normally used for disk profiles, which translates into 
Eq. (H) when rewritten to effective parameters and which 
was also used for the 2D fit. 

More stable results were reached than with the general 
exponential law, but the intrinsic properties of the r^^^ 
law still made it impossible to reach convergence in many 
cases. With n — 2 the situation improved considerably and 
using an exponential profile for the bulge, the fitting con- 
verged for all galaxies except for UGC6028. This galaxy 
has a very small bulge in a Type II profile (Freeman 1970) 
and therefore the fitting routine tends to make the bulge 
negative in order to create a hole in the disk profile. The 
fits with the exponential bulge are the most stable; con- 
vergence is reached more often and the fit results are less 
sensitive to initial values with the exponential bulge than 
with the other tested profiles. 



The tests on the artificial images in Sect. 3.4 showed 
that one of the main sources of error in the bulge and disk 
parameters was the uncertainty in the sky background 
level. The maximum errors in the fit parameters due to 
this uncertainty were calculated for each galaxy using the 
sky errors determined in Paper I. All ID fits were repeated 
with the maximum sky error estimate added to and sub- 
tracted from the Imninosity profile. These error estimates 
will be used in the next sections. 



4-. 2. Profile comparison 

The different decompositions of eight galaxies are dis- 
cussed in more detail as it will give some insight in 
the problems dealing with bulge/disk decompositions. 
Figure ^ shows some typical best and worst cases of 
four different decomposition models on B passband pro- 
files/images. Images of these galaxies can be found in Pa- 
per I. 

UGC 89 is probably the clearest example in the sample 
of a galaxy with three distinct components, which can 
only be fitted correctly using a 2D fitting technique. In 
Fig. ^ we see that in the r^/'* law fit the bulge fills up 
the bar region, but also replaces part of the disk. The 
central surface brightness of the disk is probably too 
faint and the scalelength too large in this fit. When 
we use an r^^^ or an exponential bulge, the disk tries 
to fill up the bar region, making the fit to the outer 
exponential part of the profile worse. The disk param- 
eters are obviously incorrect for these fits, but now the 
surface brightness is too faint and the scalelength too 
small. Only including a bar yields a satisfying result, 
also displayed in Fig. |l|. 

UGC 438 is a galaxy with a small bulge and a distinct 
exponential disk with some enhanced star formation 
in the spiral arms near the center. All fitting meth- 
ods seem to be equally justified, but notice the very 
extended bulge in the case of the r^/^ bulge profile. 

UGC 1455 is a symmetrical galaxy with a small bar/oval 
component in the center. All fits seem equally justifi- 
able as long as we assume that the oval component is 
part of the bulge. The /io becomes gradually fainter 
going from r^/^ to the exponential bulge. If we assume 
that the bar is part of the disk, which is likely as some 
spiral arms start at the ends of the bar, the 2D fit is 
the best fit. 

UGC 1551 has a low surface brightness disk with some 
fiocculent star formation and a very small bulge. The 
j.i/4 profile did not converge for the B passband 
(though it did for the other passbands). Even though 
convergence was reached with the other fitting func- 
tions, the bulge results can hardly be called reliable; 
the small bulge is very hard to fit. 

UGC 6077 is a typical case of a Type II profile, where 
the exponential outer profile does not continue all the 
way inward before the bulge takes over the luminosity 
in the central part of the galaxy. Morphologically these 
galaxies are in general barred with strong spiral arms 
forming at the end which are tightly wound to form a 
ring. Sometimes the bar is less pronounced. Moving to 
longer wavelengths the bump becomes less pronounced 
as the enhanced star formation doesn't show up so 
clearly at those wavelengths. None of the tested fitting 
models is able to decompose this type of galaxies well. 
One has to adjust the exponential law of the disk to 
contain an inner hole or one has to model spiral arms 
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Fig. 2. Some examples of the fits made to the B passband surface brightness profiles. The squares indicate the measured profile, 

the dashed line represents the fitted disk, the dotted line the fitted bulge and the dashed-dotted line the bar. The full line is the 
sum of the different model components. The method to extract the model luminosity profiles from the 2D model images was 
the same as used on the real data. 
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Fig. 3. The residuals A/i between data and model for the galaxies plotted in Fig. |^. The dotted line indicates the model with 
the r^^'^ law bulge, the dashed line the r^^^ law bulge fit, the dashed dotted line the exponential bulge fit and the full drawn 
line the 2D fit method. Note that generally the 2D fit method results in smaller residuals, but that the difference between the 
models is small compared to the intrinsic fluctuations in the luminosity profile. 
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Fig. 4. The average reduced values per morphological type for the ID fit with different bulge models. 



in the disk. None of the fitting techniques are correct 
(the r^/"* law fit did not even converge for this galaxy) 
and the central surface brightnesses and scalelengths 
found are just an indication of the extent and surface 
brightness of the disk. 
UGC 6445 , a galaxy with a small oval component, has 
a less pronounced Type II profile. Because of a slight 
shift in steepness of the luminosity profile in the outer 
region of the galaxy, the r^l^ fit will tend to dominate 
in the central and outer region. Consequently bulge 
and disk are equally important over the whole extent 
of the galaxy in such a fit. Inclusion of a bar or oval 
component immediately diminishes the bulge contribu- 
tion at large radii, even if one would fit an r^/* bulge 
profile. Clearly the 2D fit is the best option for this 
galaxy. 

UGC 10437 is an example of one of the galaxies in the 
sample with a very low surface brightness disk. The 
fits made to this kind of galaxy are in general very un- 
stable as the disk disappears rapidly in the sky noise. 
This means that the profile that can be fitted only ex- 
tends over a few scalelengths. The H and K passband 
observations especially could pose some problems as 
these showed hardly anything else other than the bulge 
above the sky noise (see also UGC 334, UGC 628 and 
UGC 9024). The fits with r^/* and r^/^ bulge profiles 
only converged to unphysical values for this galaxy, the 
exponential ID bulge and the 2D fit yield nearly the 
same result, as might be expected for a nearly face-on 
galaxy. 

UGC 12776 is a galaxy that consists mainly of a bar, 
with two arms originating at the ends. Similar galaxies 
are UGC 2368 and UGC 10083. Obviously a bar should 
be included in the models to fit these galaxies, but 
we have to keep in mind that in these cases the bar 
really seems to belong to the disk. The central surface 
brightness of such a galaxy is probably best described 



by the sum of the contributions of the disk and the 
bar. 

Figure ^ shows the residuals between models and data. 
This figure shows that including a bar in the fitting rou- 
tine can indeed improve the overall fit, but also shows 
that the residuals are dominated by fluctuations in the 
profiles. These fluctuations are not due to signal-to- noise 
problems, but are intrinsic to the light distribution of the 
galaxy. This makes it difflcult to decide which bulge model 
is best, as the change in reduced is small compared to 
the reduced value itself. One would need additional 
components to fit the fluctuations, but this is beyond the 
scope of this work. 

Comparison of values 

An exponential bulge was used in the 2D fltting technique, 
which was motivated by the ID test results. One such test 
is the comparison of the reduced values for the differ- 
ent bulge models. Figure ^ shows the different reduced 
values as function of morphological type. Note that the 

values are smaller in the K passband than in the B 
passband and that the differences between the x^ values 
of the different bulge models are small compared to these 
X^ values themselfs. 

For galaxies with types later than Sb, the flts with the 
exponential bulge give the smallest residuals. For galaxies 
in the range from Sa-Sb, the r^/^ law bulge flts give the 
best results. In this range, the r^/'* law and exponential 
law bulge give comparable results in reduced x^ values. 
Therefore one could propose an intermediate transition 
type bulge for early spirals between the r^^^ law ellipticals 
and the exponential bulges of late-type spirals. This could 
be analogous to the pure elliptical systems, where low lu- 
minosity systems are less centrally concentrated than high 
luminosity systems (Young & Curie 1994 ). 

The results of the study by Andredakis & Sanders 
(1994) are conflrmcd using this sample. Exponential 
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Table 4. The median of the fit parameter errors for all galaxies, where the maximum error in the sky background subtraction 
was used to calculate the uncertainty in the fit parameters. The errors in /io and /ic are in mag arcsec"'^, the relative errors in 
h and were defined by Ah/h and Arc/re. 



Band median error in /io median error in h median error in /ie median error in rc 





marking 


^1/4 


^1/2 


exp. 


marking 


^1/4 


^1/2 


exp. 


^1/4 


j.l/2 


exp. 


^1/4 


^1/2 


exp. 


B 


0.02 


0.07 


0.06 


0.05 


0.037 


0.044 


0.033 


0.029 


0.15 


0.06 


0.02 


0.50 


0.26 


0.14 


V 


0.03 


0.08 


0.07 


0.05 


0.033 


0.035 


0.030 


0.025 


0.18 


0.06 


0.02 


0.33 


0.21 


0.12 


R 


0.02 


0.07 


0.06 


0.04 


0.027 


0.034 


0.026 


0.023 


0.11 


0.05 


0.02 


0.40 


0.18 


0.09 


I 


0.03 


0.11 


0.10 


0.07 


0.042 


0.030 


0.026 


0.021 


0.18 


0.08 


0.03 


0.43 


0.15 


0.08 


H 


0.04 


0.16 


0.20 


0.12 


0.041 


0.039 


0.028 


0.023 


0.28 


0.12 


0.04 


0.43 


0.15 


0.09 


K 


0.04 


0.17 


0.14 


0.11 


0.049 


0.032 


0.025 


0.021 


0.27 


0.09 


0.04 


0.39 


0.16 


0.09 



Table 5. The mean relative change in disk parameters due to different techniques for all galaxies where the fit routine converged 
to physical values for all techniques. All mean changes are relative to the "marking the disk" method, with A^o in mag arcsec"^, 
d = 2(/imark — Mother )/ (/imark + ^othor) dimensioiilcss. Tho crrors arc standard deviations. 
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bulges are statistically at least as justified as r^/'' law 
bulges. When the galaxy image does not contain enough 
information to determine the shape of the bulge, the ex- 
ponential bulge is preferred for galaxies with classification 
later than Sb. For early-type spiral galaxies an r^/^ law 
bulge is a good alternative. For consistency reasons, only 
the results obtained with the exponential bulge will be 
used in subsequent papers when analyzing this large sam- 
ple of galaxies. 

Comparisons of errors due to sky uncertainties 

The tests on artificial images discussed in Sect. ^ indicated 
that the uncertainty in the sky background level is, once a 
bulge model has been chosen, the largest source of error. 
As explained before, the maximum errors due to incorrect 
sky background subtraction were only calculated for the 
ID methods. The median relative errors for each parame- 
ter and each method are listed in Table ^. Two trends can 
be seen in this table: 

1) the median errors in the surface brightness parameters 
increase going from the B to the K passband, while the 
uncertainties in the scale parameters decrease and 

2) the median errors in all parameters decrease moving 
from the fits with the more centrally peaked r^/^ law bulge 
profile to the fits with the shallower exponential bulge pro- 
file. 

A change in sky background level affects mostly the 
faint outer regions of the galaxy profiles. The r^/* law pro- 
file offers often a greater contribution to the outer parts 



of galaxies than the exponential bulge profile and will be 
more affected by a sky level change. The parameters re- 
sulting from the fit with the exponential bulge are the 
most stable and reproducible according to the results in 
Table ^ which is independent of the conclusion whether 
this bulge model is correct or not. 

4-.5. Comparison of the resulting disk parameters 

As a final point in the comparison of the different bulge 
models and fitting methods the resulting parameters are 
compared. Table || lists the average differences in the disk 
parameters using the different methods. The disk param- 
eters of all different models are compared to those of the 
"marking the disk" method. The change in central surface 
brightness A/ig = /J-o, mark-Mo, other is on average the largest 
for the r^/^ law bulge and small for the exponential bulge 
in both the ID and 2D case. The central surface brightness 
of the disk is lower for the r^/^ and r^/** law bulges as part 
of the disk light is replaced by the light of these more ex- 
tended bulges. The average relative change in scalelength 
{d = 2(/ii„ark - /iothcr)/(/imark + Kthcr)) is a fcw percent 
at most for all models. The scalelengths become on av- 
erage a little bit smaller. The standard deviations of the 
differences in /iq are much larger for the r^^^ and r^/^ law 
bulges than for the exponential bulges. Apparently the 
methods using exponential bulges yield more often disk 
parameters that are in accordance with the "marking the 
disk" results, which represent the human intuition. 
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Fig. 6. A comparison between the K passband parameters obtained with the ID and the 2D methods. An exponential bulge 
profile was used for both the ID and the 2D method, a) Central surface brightness of the disk, b) effective surface brightness 
of the bulge, c) scalelength of the disk and d) effective radius of the bulge. The errors in the 2D parameters are the formal fit 
errors of the fit routine, the errors in the ID parameters are the maximum errors due to the uncertainty in the sky background 
subtraction. 



A comparison the results from the "marking the disk" 
fit with the 2D fit (Fig. ||), shows that the central sur- 
face brightnesses agree for most galaxies to within 0.5 
mag arcsec"^. The few exceptions to this can be easily ex- 
plained. On the one hand, the surface brightness is fainter 
in the 2D case for UGC 2368, UGC 10083 and UGC 12754, 
as a bar is making up a large fraction of the disk and 
for some low surface brightness galaxies (UGC 628 and 
UGC 12732), because the distinction between the nearly 
absent bulge and the disk is hard to make. On the other 
hand, UGC 7901 and UGC 11628 have a gradual change in 
slope of their luminosity profile, which is interpreted dif- 
ferently by the 2D fitting program than by the human eye. 
The 2D fit central surface brightnesses are much brighter. 

Figure ^ shows a comparison between the ID and 2D 
bulge and disk parameters, both obtained with the expo- 
nential bulge profile. The errorbars on the 2D results in- 



dicate the formal fitting errors of the fitting routine, thus 
the uncertainty in the minimum found by the routine. 
The errors in the ID parameters result from recalculations 
of the fit with the sky background uncertainties added and 
subtracted. 

Two points can readily be made from these errors: 1) 
the errors due to sky uncertainties are almost always much 
larger than the formal fitting errors and 2) the data points 
which are the most deviant from the equality line have in 
general also the largest uncertainty due to the sky errors. 
The errors due to sky uncertainty are apparently also in- 
dicators of the differences between both fitting methods. 

The ID and 2D disk parameters in Fig. || agree quite 
well without any systematic deviations. The effective sur- 
face brightness of the bulge shows a small but system- 
atic deviation, the surface brightnesses determined with 
the ID method are most of the times a little fainter. The 
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Fig. 5. The difference between the central surface brightnesses 
obtained with the "marking the disk" method and the 2D fit 
method (A/io = A'™''* ~ A'o^) function of fio in the K pass- 
band. The UGC numbers of the most deviant galaxies are in- 
dicated. 



comparison of the effective radii shows the largest devia- 
tions. The 2D effective radii are smaller, especially for the 
galaxies where a bar was fitted. The profile comparisons 
of Fig. H showed that these smaller effective radii are more 
realistic. 

4-. 6. The exponential bulge versus other bulge models 

The exponential bulge profile was used in the 2D fitting 
technique instead of the more widely used r^/* law profile. 
The previous tests have shown that this is justified. The 
profile decompositions look at least as good as with other 
bulge models, which is rcfiected in the reduced results. 
Table ^ shows that the exponential bulge results are least 
effected by the uncertainty in the sky level and Tabled 
shows that the exponential bulge fit reflects the "marking 
the disk" results the most. The choice for the exponential 
bulge is obvious if one further realizes that the exponen- 
tial bulge fit, contrary to the other models, almost always 
converges without fine tuning the initial parameters. One 
could propose to use r^^^ law bulges for early-type spiral 
galaxies, but for consistency reasons I will use only the 
exponential bulge law results in subsequent papers. 

5. Error discussion 

Good error estimates of the fit parameters are important 
to determine the reliability of the derived relations in the 
next papers of this series. All sources of error and the 
typical percentage of error they cause in the different bulge 
and disk parameters are discussed. 

In the decomposition methods discussed in this arti- 
cle one is confronted with three kinds of uncertainties: 1) 
uncertainties in the component profiles, 2) the uncertain- 
ties due to the fitting method and 3) measurement errors. 



The first two items give systematic effects in the derived 
parameters, but make comparisons within such a frame- 
work still possible. The last item has no systematic effect, 
assuming measurement errors are randomly distributed. 
Each of these points will be addressed separately, indicat- 
ing the order of magnitude uncertainty in the resulting 
parameters. 

1 ) Errors due to uncertainties in component profiles 

— The best choice for the description of the bulge light 
is probably the generalized exponential (Eq. (|^)), as 
it includes both the r^/'* and the exponential law. Its 
mathematical peculiarities make it here less useful in 
combination with an exponential disk profile. There- 
fore the exponential function for the bulge is preferred, 
as it gives a first moment analysis of the part of the 
bulge that is seen above the disk luminosity. In this 
way, a surface brightness and a scale size of the central 
region of the bulge can be determined. These quan- 
tities must have some physical meaning, even though 
they show only a weak relation to the ones obtained 
with r^/^ and r^/^ law bulges. The average value of fxo 
of the disk will be 0.2 mag arcsec"^ fainter using these 
bulges and the rms differences will increase to 0.4 mag 
arcsec"^ with respect to the "marking the disk" fit. 
The scalelengths of the disks are about the same for 
all bulge models used. The uncertainty in the disk pa- 
rameters for a typical galaxy in this data set is larger 
due to bulge profile uncertainty than due to sky back- 
ground uncertainty (compare results in Table |^ and 
standard deviations in Table |5|). 

— The use of exponential light profiles for disks of spiral 
galaxies is well established, but Type II profiles intro- 
duce an uncertainty. Are we dealing with a hole in the 
light distribution (intrinsic or due to dust) or do we 
have an extra stellar component at intermediate radii 
in the form of a ring of spiral arms? Morphology and 
color considerations favor the second option. Therefore 
one has to include the inner region in the fit as this is 
also part of the disk and observe at wavelengths which 
are least affected by young stellar populations. In the 
B passband, the Type II profiles can give uncertainties 
of order 20% in the disk parameters and up to 50% in 
the bulge parameters. In the K passband the effect is 
strongly reduced to order 10% for all parameters. 

2) Uncertainties induced by different fitting methods 

— The uncertainties in the "marking the disk" method 
have been discussed by Gio vanell i ct al. ( 1994 ) and by 
Knapen & van der Kruit ( |l99l|) . They showed that 
the uncertainties using this method are of order 15%. 
It has been argued that this method gives systemati- 
cally wrong resul ts (Ko rmendy 1977, Phillipps & Dis- 



ney 1983 Davies 1990 ) which will be discussed in Pa- 
per III (de Jong |l995| ). 

The decomposition of the ID profile in a bulge and 
a disk component should in principle be better than 
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the previous method, but when performed with an 
exponential bulge profile the disk parameters hardly 
change compared to the "marking the disk" results. 
Byun ( 1992 ) showed that this method gives systemat- 
ically wrong results for inclined galaxies. 

— The 2D fitting method is intrinsically the best method 
tested here. It includes the effect of inclination and 
makes it possible to introduce a non-axisymmetric 
component in the form of a bar. With respect to the 
double exponential ID fit it has little effect on the disk 
parameters and on /ic of the bulge. The estimate of rc 
of the bulge improves. The tests on artificial images 
showed that the errors in the parameters are at most 
a few percent for bulges with larger than 1". 

— One source of error that has never been investigated in 
detail is the choice of weighting function for the fitting 
routine. I have only used weighting schemes which re- 
duce relative errors between model and data, but one 
could argue that the inner part of galaxies should get 
more weight as the signal-to-noise is better in these 
regions. Likewise one could also argue that the outer 
part of the luminosity profile should get more weight 
as many more points are sampled (this is what hap- 
pens automatically in the 2D fitting technique) . There- 
fore even a single model for bulge and disk luminosity 
distribution has an intrinsic uncertainty depending on 
where one puts the most weight for the fitting accu- 
racy. It is hard to define an acceptable range of weight- 
ing functions, but the uncertainty in the disk parame- 
ters induced in this way is estimated to be at most 5% 
for most galaxies. It will be of order 10% for the bulge 
parameters. 

3) Uncertainties caused by measurement errors 

— The dominant source of error in both bulge and disk 
parameters is the uncertainty in the sky background 
subtraction. As a second order effect b/a is also im- 
portant for the estimate of /xo- The seeing estimate 
influences the determination of the parameters of rel- 
atively small bulges. Each of these errors has at most 
order 10% effect on the disk parameters and a 20% ef- 
fect on the bulge parameters. The errors will be slightly 
larger for low surface brightness galaxies, due to the 
relatively larger contribution of the sky background 
error. 

— Another source of error, independent of the fitting 
method employed, stems from the determination of the 
zero-point of the magnitude scale. In Paper I, we cal- 
culated that this was for our sample in the range from 
0.03 to 0.12 mag arcsec"^, resulting in this contribu- 
tion being of the same order of or less than the above 
mentioned measurement errors. 

6. Conclusions 

The bulge, disk and bar parameters of the 86 spiral galax- 
ies presented in Paper I were calculated. In the different 



decomposition methods that were explored, we were con- 
fronted with three kinds of uncertainties: 1) the uncertain- 
ties in the component profiles, 2) the uncertainties due to 
the fitting method and 3) measurement errors. The first 
two items will give systematic effects in the derived param- 
eters, but comparisons within such a framework are still 
possible. The last item has no systematic effect, assuming 
measurement errors are randomly distributed. 

In conclusion, the use of the 2D fitting technique with 
exponential light profiles for both bulge and disk yields 
the most reproducible and representative component pa- 
rameters. Assuming that the 2D model with exponential 
profiles for both bulge and disk is a reasonable correct de- 
scription of the global structure of a galaxy, the estimated 
maximum errors of the structural parameters are of order 
20%. These errors are dominated by the uncertainty in the 
sky background level. The errors in the disk parameters 
will double if the model for the bulge is better represented 
by an r^/^ law profile. The errors on the bulge parameters 
will then be very uncertain. The determined parameters 
will be used in the subsequent articles in these series. 

The two dimensional decomposition technique ex- 
ploited here has many advantages above other often used 
decomposition techniques. Still we are far away from full 
scale 3D models which link complete self-consistent dy- 
namics of the different components to a 2D morphological 
projection. 
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